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Abstract. In this paper we establish a necessary condition for the application of stellar population synthesis models to observed 
star clusters. Such a condition is expressed by the requirement that the total luminosity of the cluster modeled be larger than the 
contribution of the most luminous star included in the assumed isochrones, which is referred to as the Lowest Luminosity Limit 
(LLL). This limit is independent of the assumptions on the IMF and almost independent of the star formation history. We have 
obt ained the Lowest Lu minosity Limit for a wide range of ages (5 Myr to 20 Gyr) and metallicities (Z=0 to Z=0.019) from 
the lGirardi et alj l2002t) isochrones. Using the results of evolutionary synthesis models, we have also obtained the minimal 
cluster mass associated with the LLL, M""", which is the mass value below which the observed colors are severely biased 
with respect to the predictions of synthesis models. We explore the relationship between AT"" and the statistical properties 
of clusters, showing that the magnitudes of clusters with mass equal to AT"" have a relative dispersion of 32% at least (i.e., 
0.35 mag) in all the photometric bands considered; analogously, the magnitudes of clusters with mass larger than 10 x AT"" 
have a relative dispersion of 10% at least. The dispersion is comparatively larger in the near infrared bands: in particular, AT"" 
takes values between 10 4 and 10 5 M Q for the K band, implying that severe sampling effects may affect the infrared emission of 
many observed stellar clusters. As an example of an application to observations, we show that in surveys that reach the Lowest 
Luminosity Limit the color distributions will be skewed toward the color with the smallest number of effective sources, which 
is usually the red, and that the skewness is a signature of the cluster mass distribution in the survey. We also apply our results 
to a sample of Globular Clusters, showing that they seem to be affected by sampling effects, a circumstance that could explain, 
at least partially, the bias of the observed colors with respect to the predictions of synthesis models. Finally, we extensively 
discuss the advantages and the drawbacks of our method: it is, on the one hand, a very simple criterion for the detection of 
severe sampling problems that bypasses the need for sophisticated statistical tools; on the other hand, it is not very sensitive, 
and it does not identify all the objects in which sampling effects are important and a statistical analysis is required. As such, it 
defines a condition necessary but not sufficient for the application of synthesis models to observed clusters. 
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1. Introduction and motivation 

The comparison between theoretical models and observations 
is the basic procedure that allows the evaluation of our current 
knowledge of Nature; underlying this statement is the assump- 
tion that better observational data, by setting more stringent 
constraints, make such comparison more meaningful. Although 
this is almost always the case, there is a situation in which, 
paradoxically, the opposite is true, and the more the data quality 
improves, the more biased the theoretical inferences turn out to 
be. This is indeed the case of the analysis of the integrated light 
of stellar populations. The progressively higher sensitivity of 
modern instruments provides access to data from increasingly 



fainter stellar clusters, up to a point where we must begin to 
take into account the limitations posed by the discreteness of 
the number of stars in a system: due to incomplete sampling, 
both small clusters and small samples of large clusters have 
stellar mass distributions that may differ substantially from the 
one predicted by the underlying initial mass function (IMF). 
Nevertheless, most theoretical models assume that the IMF of 
clusters is completely populated, i.e., that the distribution of 
stellar masses is continuous and that all the evolutionary stages 
are well sampled. Obviously, any model assuming a continu- 
ous IMF (hereinafter, analytical models) will be correct only 
under the asymptotic assumption of an infinite number of stars. 
Hence, the validity of the comparison between the predictions 
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of synthesis models and real systems, where the IMF is not 
perfectly sampled, depends on the size of the system. 

Several works have been written that deal directly or 
indirectl y with the subjec t of sam p ling effects. F or ex - 
ample. iBarbaro & Bertelhl dl977l). IChiosi et all |l988), 
iGirardi & Bical dl993h . and iGirardi et alJ dl995l) . who 
show the relevan ce of sampling effects for the study of 
LMC clusters; ISantos & Frog ell lll997l) . who determine 
how sampling effects affect integrated near-in frared colors; 
ICervino. Luridiana . & Castanderl J2000l) and ICervino et alJ 
(2000J), who study the effects of sam pling in some observ- 
able s of young star-forming regions: lCervifio et al.l (|2p01a) 
and ICervi no & Molll d2002l) . who estimate the effects of 
sampling in stellar yields and chemical evolutionary models. 
Sampling effe cts also underlie the stud y of surface brightness 
fluctuations dTonrv & Schneider! Il988l) . a primary distance 
indicator that is based on the analysis of the variations with 
distance of the amount of stars sampled by CCD pixels: for 
example, ICantiello et al 1 d2003l) show that surface brightness 
fluctuations may suffer from a bias that depends on the density 
of stars in the image pixels. 

In almost all the preceding works, sampling effects have 
been evaluat ed by the us e of M onte Carlo simulations. 
Alternatively, Cervino et al. (2002b) proposed a formalism, 
based on the original formulation by iBuzzonil dl989l) . where 
the mean value, the dispersion, and the correlation coefficient 
of different observables are obtained analytically using a con- 
tinuously distributed IMF. The method is applied to young star- 
forming regions (f < 20 Myr) and compared with the results of 
Monte Carlo simulations, showing that in both cases the results 
are quite similar, except for colors and equivalent widths in 
clusters with a low number of stars . The m ethod is extensively 
tested in lCervifio & Valls-Ga baud (2003) for clusters with a 
number of stars between 1 and 10 3 , where it is shown that it 
reproduces the average value and the dispersion of quantities 
obtained from Monte Carlo simulations, i.e., the luminosities 
of Monte Carlo simulations are distributed around the mean 
value obtained from the analytical model, if the quantity scales 
linearly with the amount of stars in the system. However, when 
the modeled properties are logarithmic quantities or ratios, as 
in the case of equivalent widths and colors, the mean values of 
Monte Carlo simulations are biased with respect to the results 
of the analytical modeling; the smaller the system, the more se- 
vere the bias. Unfortunately the authors are not able to quantify 
this bias in an analytical way for very small systems. 



The su bject of sampling h as als o 
by, e.g.. lLancon & Mquhcinel d2000l). 
Cervifio et all d2001bl) . iBruzuall (120 02). 



been a ddressed 
Idirardil (12000), 
iGirardil J2002), 



Cervifio et all ll2002al) ICervifio & Luridianal d2002ajr| . 
Cervifio & Valls-Gabaudl d2002l) . and ICervifiol d2003alblJdh . 



However, since all these works are published in conference 
proceedi ngs, their co nsequences are not extensively explored. 
Lan con & Mouhcinel (|2000) evaluate sampling effects on 
monochromatic luminosities at solar metallicity without the 
use of Monte Carlo simulations, and quote some limits for the 
minimal initial clusters masses ensuring a relat ive erro r lower 
than 10% for some ages and luminosities. Bruzual (2002) 



(see also iBruzual & Charlotl I2003I) 1 presents Monte Carlo 
simulations in which the stochastic effects on U — B, B — V, 
V - K, and K for the LMC metallicity are presented as a func- 
tion of the initial mass of the cluster. His figures show clearly 
that there is a bias in the results of Monte Carlo simulations 
with respect to the results of analytical synthesis models; 
however, this result is not me ntioned in the t ext, possibly due 
to the limited space. Finally, Girardi (2002) presents Monte 
Carlo simulations where the effect of a continuous distribution 
in the initial cluster masses is studied. His results are more 
appropriate for the comparison with surveys of real clusters 
than those that do not consider a distribution of masses. 

In all the preceding papers, the evaluation of sampling ef- 
fects requires making assumptions on the IMF and the star for- 
mation history, a fact that limits the practical application of 
the results to real observations. In this paper we propose in- 
stead a method entirely based on observable quantities, there- 
fore independent of the IMF and almost independent of the star 
formation history, to estimate whether the colors predicted by 
synthesis models are biased with respect to real observations, 
and to establish when analytical synthesis models cannot be ap- 
plied and a statistical formulation is required. We also describe 
the relationship of this method to more sophisticated statisti- 
cal analyses. Furthermore, we suggest examples of possible ap- 
plications of the method to the analysis of observational data. 
Finally, we discuss the qualities and drawbacks of our method 
in comparison to alternative ones. 

The structure of the paper is the following: in Sect. 2 we 
define the method, and provide a quantitative evaluation of the 
quantities involved for a wide range of observational cases; in 
Sect. 3 we translate the preceding results in terms of the clus- 
ter masses, proposing a mass criterion to exclude low-statistics 
clusters from the analysis performed with synthesis models; in 
Sect. 4 we show the observational implications of this work; in 
Sect. 5 we discuss the limitations of our results; and in Sect. 6 
we draw our conclusions. 

2. The Lowest Luminosity Limit 

One of the most basic limits to the application of evolutionary 
synthesis models can be expressed by the following statement: 

The total luminosity of the cluster modeled must be larger 
than the individual contribution of any of the stars included in 
the model 

This obvious statement defines a natural theoretical limit 
that has not always been considered when models are applied 
to real observations. W hereas in the work by Tinsley (see, e.g., 
iTinslev & Gunnlfl976l) it was not necessary to take this limit 
into account, due to the observational limitations at that epoch, 
the increasing sensitivity of current instruments has reached a 
level where this limitation plays a fundamental role in the in- 
terpretation of data. 



1 A version of Bruzual & Chariot 1 2003) synthesis models can 
be obtained at http://www.cida.ve/~bruzual/bc2003 or at 
http : //www2 . iap . fr/users/charlot/bc2003 
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Defining as the Lowest Luminosity Limit (hereinafter LLL) 
the luminosity of the brightest individual star included in the 
model, we can establish a simple luminosity criterion for the 
application of synthesis models to the interpretation of ob- 
served clusters, imposing the condition that the cluster modeled 
be more luminous than the LLL. While clusters brighter than 
this limit may either be well-sampled or not, clusters fainter 
than this limit are certainly misrepresented by synthesis mod- 
els. 

It is possible to establish a luminosity limit following 
different criteria, yielding either weaker or more stringent 
constraints; but all the alternative definitions we could think 
of turned out to either lack physical meaning, or imply 
circumstances that do not occur in the modeling practice. 
Nevertheless, the degree of arbitrariness of our definition is a 
hairy problem, and we will discuss it thoroughly in Sect. [5] 

With the definition given above, the LLL is only defined 
by the isochrone used and the band under consideration; how- 
ever, its exact value at a given age is also weakly dependent 
on the star formation history. In the following we present the 
values of the LLL computed for a wide range of parameters. 
We have used the isochrones and the integrated magnitudes 
of simple stellar population models (hereinafter SSP models, 
i.e., model s that assume an in stantaneous burst of star for- 
mation) by Girardi et al.l d2002l) 2 . We consider seven different 
metallicities: Z=0.019 (solar), Z=0.008, Z=0.004, Z=0.001, 
Z=0.0004, Z=0.0001, and Z=0.0. The Z=0.0 isochrone s cor- 
respond to metal-free models by Marigo et aO feOOlll . The 
other isochrones correspond to the basic set presented in the 
web serve r cited in the footn ote, which combines the re- 
sults from Idirardi et alJ feOOOt and Idirardil fcOOlh for low- 
arid in ter mediate-mass stars, with the results by Bertelli et al. 
(1994) and lGirardi et al .(1996) for high-mass stars, and that in- 
cludes overshooting and a simplified Thermal Pulse AGB (TP- 
AGB) evolution. Additionally, f or Z=0.019, 0.008 , and .004 
we have used the isochrones by Marigo & Girardi J200l[) that 
include a more detailed TP- AGB evoluti on. Most of the atm o- 
sphere models are taken from ATLAS9 JCastelli et alJll997l) 3 . 
A more detailed description of the isochrones, the atmosphere 
model s, and the SSP models can be found in iGirardi et alJ 
(2002) and in their web server. 

To obtain the LLL we have searched for the bolometric 
luminosity of the most luminous star in any given isochrone 
LTf(t). However, since the tabulated data also provide the 
magnitudes at different bands, we have also obtained the mag- 
nitude of the most luminous star in the Johnson-Cousins-Glass 
filters U, B, V, R, I, J, H, and K: M'l} in (t), M™ n (i). 

In general, the most luminous star is also the most evolved, 
but this relation does not strictly hold for all the bands nor 
all the ages. This fact is well illustrated in Fig. ^ where the 
LLLs obtained from isochrones with a simplified TP-AGB 
treatment (right panels) are compared with the LLL obtained 
from isochrones with a detailed TP-AGB evolution (left panels) 
for models with Z=0.019, 0.008, and 0.004. Whereas M™ n B VRJ 

2 Isochrones and simple stellar population results are available at 
http : //pleiadi . pd. as tro. it/| 

3 NOVER models at http : / /cf aku5 . harvard . edu/grids . html 



depend on the TP-AGB treatment at some ages, M^"" JHK are 
almost unaffected by it. 

Figures [2 and [5] show the values of the LLL in magni- 
tudes for different ages and metallicities. The figures show that 
M m "'(t) evolves with time toward less luminous values for all 
the bands, and that blue bands become fainter more quickly 
than red bands, as expected due to the cluster evolution: the 
population becomes redder and fainter as the cluster evolves. 
These two trends combined imply that, as the cluster ages, 
the LLL tend to decrease (i.e., become less stringent) in all 
bands, doing so more rapidly at blue wavelengths than at red 
wavelengths. It is also interesting to compare the evolution of 
M'™(f) with the various M min (i) at different bands. It can be 
seen that the larger the metallicity, the more similar is the evo- 
lution of M™"(f) to the evolution of M mm (t) in red bands, as a 
consequence of the increasing fraction of bolometric light that 
goes into red filters. The bottom panel in Fig.[2]compares the 
LLL values for V and K and two extreme metallicities. It can 
be seen that M"™ is almost metallicity-independent. 

These figures show the LLL values for the case of single 
isochrones, that correspond to SSP models. However, the LLL 
can be easily obtained for different star formation histories. For 
example, let us consider a two-burst system with ages ti and 
ti, associated to two minimum magnitude values M m '"(t\) and 
M mm (ti): then the LLL of this system is simply the minimum 
of M min (t\) and M m, "(f 2 ). In the case of a cluster with con- 
stant star formation history and age f, the LLL at each band 
is given by the maximum of the luminosity of SSP models in 
the age range between and f, which, for the bolometric lumi- 
nosity, can be roughly approximated by the ZAMS luminosity 
of the most massive star included (this result is not exact be- 
cause the bolometric luminosity of a star briefly increases af- 
ter the ZAMS). At other bands the maximum luminosity can 
be reached much later, thus no easy recipe can be given. As a 
final remark, note that some evolution in the metallicity is ex- 
pected for any star formation history different from an instanta- 
neous burst, and this effect should be included in the modeling 
in a self-consistent way (see Schulz et al. 2002, and references 
therein as an example) 4 . 

Summarizing, the most luminous star in one band is not 
necessarily the most luminous star in the rest of the bands, 
neither is it the star with the largest bolometric luminosity. 
As a consequence, the LLL depends not only on the age and 
the metallicity (i.e., the assumed isochrones and atmosphere li- 
braries), but also on the band considered. The dependence on 
age and metallicity is weaker for the near infrared bands than 
for the optical bands. In particular, the LLL for K is almost 
metallicity-independent. 



3. Minimal initial cluster masses 

In this section we will translate the concept of LLL into an 
equivalent formulation in terms of mass. 



4 The reference corresponds to the synthesis 
code galev and the model results are available at 
http : / /alpha . uni - sw . gwdg . de/ ~galev/ 
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Let us recall that evolutionary synthesis models are based 
on the convolution of isochrones with the IMF and the star for- 
mation history. For the case of a SSP, the mean luminosity in 
a given band and at a given age, l ssp (t), results from the sum 
of the luminosities of individual stars at the corresponding age 
as given by the isochrone, /,(m,, f), weighted by the number of 
stars with initial mass m, as given by the IMF, w^m,). If the sum 
of the Wj values is normalized, as usual, to 1 M Q transformed 
into stars from the onset of the burst, the resulting luminos- 
ity will also be normalized 5 . The total luminosity of a modeled 
cluster, L c,us (t), is directly proportional to the initial mass trans- 
formed into stars in the cluster, At: 



L dus (t) = M x Wiimd kim, t) = M X r sp (t). 



(1) 



Then, for a given age and metallicity, we can obtain the 
total initial mass transformed into stars from the observed lu- 
minosity L clus (or M c,us expressed in magnitudes) and the cor- 
responding normalized value of l ssp (t) (or m ssp ): 



clus 



M = 



[ssp 



2.5 log M = m s 



■M c ' 



(2) 



where we have dropped the explicit reference to t to simplify 
the notation. Now, imposing that M clus - M mm (t), we can ob- 
tain the initial cluster mass for each band and age, M m "'(t), for 
which the total luminosity of the cluster simulated by a SSP 
model equals the luminosity of the most luminous star in the 
band, M min (t): 



2.5 log M' m "(t) = m ssp (t) - M mm (t) 



(3) 



The superindex min reminds that below this limit the clus- 
ter cannot be modeled by means of a synthesis model. Note 
that Af" ! "(f) depends on the age and the band, but also on the 
IMF and the star formation history, since integrated quantities 
depend on them. Here, we have used the int egrated magnitudes 
of the SSP models bvlGirardi et all (J2002), which assume the 
IMF bv lKroup a (2001) in its corrected version (his Eq. 6), and 
a total SSP initial mass equal to 1 M Q in the mass range 0.01 - 
120 M . 

The results for different ages and metallicities are shown in 
Fig. |3] which shows that M™" l]HK are almost metallicity in- 
dependent except during the first stages of the evolution of the 
cluster. In the case of optical bands, the lower the metallicity, 
the larger M™ l B VR r The value of M""" spans three orders of 
magnitude depending on the band, and it takes a value as large 
as 10 4 M or even more for the case of near infrared colors. A 
first conclusion we can draw from these figures is that these val- 
ues are so high that many observed clusters certainly fall below 
this limit, and their properties can by no means be reproduced 
by synthesis models. This fact is too often overlooked, and we 
will try to emphasize it repeatedly throughout the paper. 



3 In the following we will use the lower case and/or the super-index 
ssp to refer to normalized quantities obtained by SSP models, and the 
upper case for absolute (denormalized) quantities. 




C02 t# l0% ; t < 10 Myr) 
LM C*i 0% ) , 
W94 (M. 10% ; t > 1 Gyr) 



GOO (LLL, Solpeter) 
G02 (LLL, Kroupa) 



Age [yr] 

Fig. 4. M mm for V (solid line with triangles) and K (dashed 
line with diamonds) assuming solar metallicity a nd diff erent 
IMF and SS P models : SSP m odels bv lGirardi etaTI d2002t) with 
the I MF bvlKrouoal J2001I) (G02: bold line), and SSP mod- 
els by Girard i et alTll2000l) with Salpeter's IMF (GOO: narrow 
line). The figure also sho ws the minimal cluste rs masses that 
ensure o~ L /L = 10% from lCervino et alJ J2002bl) models (C02: 
big open symbols until age of 10 Myr), Worthevl (119941) mod- 
els (W94: big open symbol s with ages larger than 1 Gyr) and 
lLancon & Mouhcinel (|2000i) models (LM00: small open sym- 
bols). 



In order to examine the influence of the IMF on this re- 
sult, we have also obtained the LLL from the isochrones by 
iGirardi et alJ d2000h- which ha ve been computed assuming 
a Salpeter IMF (Salpeter 1955) in the mass range 0.039 to 
100 M n from 63 Myr to 17.8 Gyr. For a comparison with 
iGirardiet alJ |2002) we have renormalized these values to the 
mass range 0.094 - 120 M in such a way that the fraction 
of mass in the range 1 - 120 M for the Rroupa's and the 
Salpeter's IMF is the same. The comparison between the M m,n 
values obtained from the two sets of isochrones and SSP mod- 
els is shown in Fig.0] With this normalization, the M mm values 
obtained using the models by Gira rdi et al, (2000) coincide for 
V with those obtained using the models by Girardi et al. ( 2002) 
at ages smaller than 80 Myr, and for K in all the age range 
in common. The differences in V for t > 80 Myr are due to 
small differences in the LLLs computed with the different set 
of isochrones. 



3.1. Minimal initial cluster mass and the estimation of 
sampling effects (theoretical point of view) 

Neither the computation of LLL nor that of Al™'" allow, by 
themselves, an evaluation of the extent of sampling effects in 
observed clusters: they just provide an easy-to-use cutoff crite- 
rion to discriminate clusters with severe sampling effects from 
clusters with more moderate or negligible sampling effects, 
without providing a quantitative tool to estimate their statistical 
properties. 
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In spite of this limitation, At""" is intrinsically related to 
sampling effects, and it is worth studying the relation between 
the information provided by At""" and the statistical properties 
of clusters. In such a way, it will be possible to obtain hints on 
the necessity to account for sampling effects in the interpreta- 
tion of observed data, even without a proper statistical formal- 
ism. 

To study the relation between the At""" value and the sam- 
pling effects, we have derived the values of the cluster masses 
associated to a relative dispersion of 10% in the luminosity of 
a band, Atio%; note that a relative dispersion of 10% in the lu- 
minosity means, at zero order approximation, cr — 0. 1 mag in 
magnitude 6 . 

To cover a wide age range, we have used the results of three 
different synthesis codes. Fo r young ages, we have used the 
solar metallicity results from ICervino et all J2002bl) models 7 . 
These models have been computed using the Geneva evolution- 
ary tracks with a Salpeter IMF in the mass range 2 - 120 M . 
In these models, sampling effects are evaluated by means of an 
effective number of stars that contribute to a given observable, 
N ssp , which is defined by: 



1 cr 2 (l ssp ) XWilf 



(I^p) 1 Q> ; /,-) 2 ' 



(4) 



as first derived bv lBuzzonil dl989l) 8 . Since the tabulated N ssp 
values are normalized to At, the value of Atio% is trivially ob- 
tained imposing that <x L /L = 0.1; therefore, the absolute effec- 
tive number of stars associated to Atio% is N — M ssp xM\o% = 
100. 

For intermedi ate ages, w e have used the results quoted by 
lLancon & Mouhcinel (|2p00) that assume solar metallicity and 
a Salpeter IMF in the mass range 0.1 - 120 M . The mod- 
els ha ve been computed with the p opulation synthesis code pe- 
gase 9 JFioc & Rocca-Volmerangell997h. 

For old stellar populations (f > 1.5 Gy r ) we h ave used the 
solar metallicity SSP models by Worthey (1994) 10 computed 
for a Salpeter IMF in the mass range 0.1-2 M Q and normal- 
ized to At = 10 6 M . This author does not compute directly 
sampling effects, but they can be inferred from the fluctuation 
luminosities, l ssp . These l ssp (or m ssp expressed in magnitudes) 
are defined as: 



r p = 



(5) 



and they are computed for the evaluation of surface brightness 
fluctuations. From Eqs.E|and|5]it is found that M ss Pxl ssp = l ssp 

6 A more detailed analysis with Taylor expansions to second or- 
der shows a small bias of 0.005 mag and an unbiased cr = 0.11. See 



ICervino & Valls-Gabaud 120031) for more details. 

7 The models used here do not include the nebular contribu- 
tion to the photometric bands. The complete set of models for 0.1 
to 20 Myr, including also the nebular contribution, is available at 
http : //www . laef f . esa . es/users/mcs/SED 

8 The outputs of A. Buzzoni synthesis code are available at 
http : / /www . merate . mi . astro . it/~eps/home . html 

9 Available at http : / /www2 . iap . fr/user s/f ioc/PEGASE . html 

10 Available at http : //astro .wsu. edu/wor they/ 



(see lBuzzonil 1 9931 for a general description of the relation be- 
tween N ssp and the brightness fluctuations). The corresponding 
N ssp can be obtained as a function of l ssp and l ssp (or m" p and 
m ssp in magnitudes) using the following formulae: 



N ssp x 10 b 
2.5 log N + 15 



l ssp 

Jssp ' 



(6) 



where the factor 10 6 (or 15 = 2.5 log 10 6 in magnitudes), is 
used to renormalize his tabulated data to Al = 1 M Q following 
a Salpeter IM F slope in the mass range 0.1 - 2 M B . 

Note that l Lancon & Moun^ <& and I Cervino et alj 
(2002b) quote the dispersion for the monochromatic luminosi- 
ties at the effective wavelength of the band (see King 1952, for 
the definition of effective wavelength) whereas Worthev ( 1994) 
(Eq. |5} gives the dispersion of the integrated luminosity of the 
band. In spite of this difference, these results can be directly 
compared: in fact, wide-band luminosities can be estimated, 
with an accuracy of 2 - 3%, by multiplying the monochro- 
matic luminosities at the effective wavelength by a constant 
value, which represents th e absolute flux density in the band 
(lKin dh 952t I.Tohnsonl 1 96(h : since this transformation only im- 
plies the multiplication by a constant, which cancels out in 
the computation of the variance, the relative dispersion of the 
monochromatic luminosity is equal to the relative dispersion of 
the luminosities in the band obtained from the exact integration 
of the spectrum over the filter resp onse. Therefore, the re sults 
bv lLancon & Mouhcinel b.QQ(h and lCervifio eTaD d2002bh and 
those by Worthev ( 1994) quoted above can be directly com- 
pared. 

In all the cases we have renormalized the resulting values 
to a Salpeter IMF in the mass range 0.094 - 120 M . The re- 
sults are shown with open symbols in Fig. |4] A first compar- 
ison among the M\o% results show that they are quite con- 
sistent with each other, with some differences that can be at- 
tributed to the differen ce between the N for malism and the 
method applied bv lLancon & Mouhcinel(l2000l) . The figure also 
shows that the evolution of Al""" is quite similar to the evo- 
lution of Atio%. The differences in the two curves range be- 
tween 0.98 and 1.5 dex (i.e., factors between 8 and 30 in initial 
cluster masses) depending on the age and the band. Therefore, 
for the assumed IMF and rounding up numbers, we can say 
that At""" is always at least a factor 10 below Atio%; this im- 
plies an N value lower than 10 and relative dispersions larger 
than 32% {cr > 0.35 mag). For these values of N, Poisson 
statistics produce non negligible probabilities of zero effective 
so urces, and hence the presence o f biases in colors, as shown 
in lCervino & Valls-Gabaudl(l2003l) . 

The relation between N and the occurrence of dispersion 
and of a bias can be easily understood in the following terms. 
Let us take as an example the case of V — K, and assume 
that L^ iL corresponds to the luminosity of a star in the Red 
Supergiant (RSG) phase. Let us also assume the case of a 10 
Myr old burst with solar metallicity, where, according to ana- 
lytical models, more than 90% of the luminosity in K is due 
to RSGs. For this case, comparing the corresponding N ssp {Lk) 
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value with the (normalized) number of RS Gs, n s J? r , i t is found 
that n^ c ~ 0.9 x N ssp {L K ) (see lCervifio et al.ll2002bl and the 
web server mentioned above). For simplicity we will assume 
that N(L K ) - Nrsg-> the absolute number of RSGs. Finally 
let us assume that N follows Poisson statistics, or, in terms 
of RSGs, that the number of such stars in different clusters is 
distributed following a Poisson distribution with a mean value 

Nrsg- 

Since the contribution of RSG stars has a small influence on 
V and a large influence on K, clusters in a mass range in which 
variations of +1 in the number of such stars are relevant will 
have colors considerably redder (dominated by an excess of 
RSGs) or bluer (due to a deficit of RSGs) than those predicted 
by analytical synthesis models. Furthermore, if the mass of the 
cluster is such that 1 < N < 10, according to Poisson statistics 
there is a fair probability that the cluster has no RSGs at all. 
In this last case (V - K) clus will be more similar to the colors 
of main sequence (MS) stars than to the colors of SSP models 
(i.e., there will be an excess of blue clusters in a survey of clus- 
ters with this mass, age, and metallicity). For these values of 
N the dispersion in the colors will be the largest. Finally, if N 
takes values lower than 1, there will be an important fraction 
(or, in extreme cases, even a majority) of clusters without RSG 
stars, and then, the mean value of the observed color (V-K) clus 
will be biased with respect to the resulting color (V - K) ssp of 
a synthesis model. On the other hand, the dispersion will de- 
crease since the range of possible (V—K) values will be smaller 
with respect to the case of larger N values. 

The situation for the case of N — 5.5 is illustrated in Fig. 
[5] which shows the results of 10 4 Monte Carlo simulations for 
clusters with 10 3 stars in the mass range 2-120 M . The ana- 
lytical value of (V - K) dus is shown by a vertical dashed line, 
and its distribution in the set of simulated clusters is shown by 
the bold solid line. Note that the mean cluster mass of these 
simulations is M = 2 x 10 4 M , a value larger than the mini- 
mal mass value M""" = 1.5 x 10 4 M : this fact, and the con- 
siderable width of the (V - K) clus distribution, confirms that 
a considerable dispersion in the observables is still expected 
for cluster masses larger than Af""\ up to values as large as 
10 x Af™ (Figure^. According to Poisson statistics, there is a 
0.4% probability of finding a cluster without any RSGs. Indeed, 
there is a small accumulation of simulations with (V - K) clws 
values around -0.6 mag, and the number of such clusters is 
about 40 (i.e., 0.4% of the total). Note also that the distribu- 
tion is negatively skewed, i.e., it tends to cut off sharply in the 
red and extend toward the blue. This means that either there 
is a deficit of RSGs, or that the RSG stars in the cluster have 
luminosities lower than the one that defines the LLL. 

This statistical interpretation depends in fact only on the 
N value, independently of it being related to a physical num- 
ber of stars (see the figures presented in lBruzualll2002l as an 
example). In general, the distribution of colors will be skewed 
toward the band with larger N ssp , i.e., toward the blue in the 
present example. 

Note that this interpretation has been done in terms of clus- 
ter masses (or, equivalently, N): the mass and the age of the 
cluster are fixed, and the dispersion in the luminosities is pro- 
duced by the random differences in the stellar mass spectrum. 



However, the cluster mass is not an observable. So a different 
approach is needed to deal with the observational problem, and 
it is more useful to use directly the LLL, as we show in the 
next section, where the remaining features of Figure|5]will be 
discussed. 

4. Applications of the Lowest Luminosity Limit 
(observational point of view) 

Up to this point, we have discussed the statistical properties of 
real clusters, trying to answer the following theoretical ques- 
tion: What is the statistical dispersion to be expected in the ob- 
servables of clusters with given mass (or number of stars) and 
age? Although answering this question surely provides a deep 
insight in the analysis of stellar clusters, it must be kept in mind 
that the observational approach to this problem differs from the 
theoretical vision: a simple reason is that, when observations 
are made, neither M nor the age of the observed clusters are 
known. The observational question can be instead put in these 
terms: Given an observed value of the luminosity, which are 
the distributions of At age, and metallicity consistent with the 
observations? 

Unfortunately, as we have repeatedly stated earlier in the 
paper, the LLL method is not a sophisticated tool, and its reach 
is very limited. More in general, the last question cannot be 
addressed with the current theory available, and a more elabo- 
rated theoretical study of this subject is needed; to this respect 
we want to remind again the work from Girardi (2002) as the 
most plausible direction to which the theoretical evaluation of 
the dispersion must be focused. However, the concept of LLL 
is powerful enough to allow us some applications to real ob- 
servational problems. Since one of the constraints imposed by 
observations is the existence of a luminosity limit, in this sec- 
tion we will explore the consequences of having a cutoff in 
luminosity in surveys of star clusters. 

In Figure [5] we described the V — K distribution of a sam- 
ple of model clusters with fixed number of stars, and noted that 
such distribution is negatively skewed. Note that a fixed num- 
ber of stars correspond roughly to a fixed cluster mass. Now, 
let us consider only the subset of clusters with a luminosity 
larger than the LLL in K, £^ iL . This subset is shown by the 
light-shaded histogram, and it can be seen that its mean color 
is somewhat redder than the predictions of SSP models. The 
key point here is that this behavior indicates that the statistics 
of the clusters is low: in fact, at t - 10 Myr an important frac- 
tion of the luminosity in V is provided by MS stars, whereas 
the luminosity in K is completely dominated by RSGs, which 
are intrinsically scarcer than MS stars. Therefore the spread in 
K among the clusters of the sample is comparatively large, and 
a cutoff in luminosity, leaving out the faintest clusters, sensibly 
alters the mean K value of the remaining subset. Schematically, 
we can say that in L K limited samples the mean K luminosity 
increases and the V luminosity is barely affected, thus V - K 
increases. This is confirmed by the application of the more re- 
strictive constraint Lf^ s > 2 x L L ^ L ", which is shown by the 
dark-shaded histogram: the number of clusters that fulfill this 
constraint is even lower, and they are redder than the predic- 
tions of SSP models. 
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Fig. 5. Probability density distribution for V-K obtained from 
10 4 Monte Carlo simulations of 5.5 Myr, with 10 3 stars in the 
mass range 2-120 M (which corresponds to vV(=2xl0 4 in the 
mass range 0.094-120 M ). The V-K distribution for differ- 
ent cutoffs in Lk is shown by different shadings (see text). The 
distribution for clusters with 10 2 stars L K > L^ Li is also shown. 



So far for clusters with (roughly) the same mass. Now, let 
us consider the behavior of simulated clusters with different 
masses. To this aim, we will assume in the following that the 
distribution of cluster masses follow s the law if/(M) oc /VL 2 , 
as proposed by Zhang & Fall ( 1999). Whether this mass dis- 
tribution correctly represents real clusters does not concern us 
for the moment: for the sake of the argument, it is enough to 
assume just any law. In fact, at the end of this section we will 
mention a possible caveat of this particular law, which might 
possibly imply a bias in the method used to derive this law. We 
have performed 10 4 Monte Carlo simulations of clusters with 
10 2 stars (/VI = 2 x 10 3 M ), which have a mean mass roughly 
1/10 that of the simulated clusters discussed up to now. Given 
the mass distribution law assumed, we have multiplied each bin 
by 100, to reproduce the expected number of clusters with this 
mass compared with clusters with mass ten times larger. 

The dotted histogram in Figure |5] shows the distribution of 
those of such clusters that also fulfill the condition L^" 5 > L^ Li . 
Note that the histogram peaks around (V - K) clus = 4.2, i.e. 
these clusters have extremely red colors; these behavior is con- 
sistent with the interpretation of being an effect of sampling, 
which in clusters of 10 2 stars is much more severe than in clus- 
ters of 10 3 stars. 

This argument can be repeated for cluster samples of any 
mass value. If we extrapolate the result to a continuous clus- 
ter mass distribution, we readily realize that the resulting his- 
togram will have a cutoff in the blue (small V-K values), 
and a long tail in the red (large V-K values): i.e., it will 
be positively skewed, contrarily to the histogram of a complete 
(non-luminosity limited) sample of clusters with a fixed mass. 
Generalizing to different colors, the observed distribution of 
colors in a luminosity-limited sample will be skewed toward 
the band with the lowest N ssp value. 



Two main conclusions can be drawn from these results. 
First, the shape of the color distribution of a luminosity-limited 
sample of clusters may be used to constrain the underlying 
cluster mass distribution, if other cluster parameters, such as 
age and metallicity, are known. To this respect note the par- 
ticular shape emerging from an extrapolation of our exam- 
ple to more mass values is just a consequence of having as- 
sumed a particular mass distribution law: in an observed sam- 
ple, the shape may in principle be different. Second, these re- 
sults also show that the LLL is an useful (probably, too con- 
servative) criterion to detect the clusters with severe sampling 
effects. Indeed, in our examples the color of the subsets with 
L c £" s > L^ LL is different from that of the complete set, indicat- 
ing that the sample suffers from sampling effects: had the mean 
mass of the clusters in the sample been sufficiently larger, no 
clusters would have been excluded by the luminosity criterion. 
Note also that the fact that the application of the luminosity 
cutoff changes the mean color of the sample implies that all the 
clusters of the sample suffer from incomplete sampling, and not 
only those that are excluded: that is, having a luminosity larger 
than the LLL is a condition necessary but not sufficient for the 
meaningful application of synthesis models; Equivalently, the 
LLL criterion detects some, but not all, of the clusters with poor 
statistics. 

To conclude this section, a note of caution about the deter- 
mination of the mass distribution law by Zhang & Fall ( 1999): 
these authors derived their law by considering only clusters 
with My brighter than -9 to avoid contamination of the sample 
by individual stars, hence applying a selection criterion that is 
more or less the observational counterpart of the LLL method. 
However, the My'" we obtain here for the age range they con- 
sider (2.5 < t < 6.3 Myr) lies between -10 mag and -9 mag: 
therefore, their analysis could possibly be affected by the sam- 
pling effects and the bias we are discussing. 



4.1. The distribution of Globular Clusters 

A clear example of a skewed distributi on is given by the sam- 
ple o f Globular Clusters (GCs) by (Gebhardt& Kissler-Patig 
1999), a fact suggesting that the properties of these objects 
might be affected by sampling effects. This statement might 
seem surprising, for GCs are the paradigm of well-populated 
objects; stars in GCs occupy the most populated part of the 
IMF, so one would naively conclude that they cannot be pos- 
sibly affected by sampling effects. However, the importance of 
sampling effects does not only depend on the absolute number 
of stars in a given mass range, but also on the evolutionary time 
scale considered: indeed, globular clusters also contain bright 
stars in low-populated evolutionary phases, such as the RGB 
and AGB phases. Hence sampling effects may also show up at 
old ages. 

Let us illustrate the prob lem by th e application of the LLL 
to GCs in NGC 5128. lReikubal J200ll) presents detailed pho- 
tometry of GCs in NGC 5128. We have used the result from 
this author because her plots can be easily reproduced with the 
data and the indications given in the paper. She compares the 
position of the GCs in color-color diagrams with the results of 
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Fig. 6. Color-color diagram for Globular Clusters in NGC 
1528. Crosses correspond to clusters with L c ^ us / L^ LL > 30 (f= 
1 Gyr); triangles to If^/Lf 1 e [10, 30]; stars to Lf ,s /L^ LL e 
[5, 10]; a nd squares to L c l " s /L^ LL < 5. Small asterisks corre- 
sponds to Worthev ( 1994) models for all the ages and metallic- 
ities quoted by the author. 



synthesis models, and notices that the clusters lie slightly to the 
right and below the model lines in the V-K vs. U—V plane. She 
mentions two possible explanations: (i) a difference between 
observations and SSP models, particularly significa nt in U - V, 
which has been hypothesized by Barmv & Huchra (2000) as a 
consequence of problems in the a tmosphere libraries used by 
synthesis models (see also B user & KurucJl978l) : (ii) an addi- 
tional offset may arise from the difference between the Bessell 
and the Johnson U -band transmission curves. These facts might 
explain the offsets between SSP models and the mean color of 
the observations, but they do not explain neither the observed 
dispersion nor the shape of the distribution. Let us now study 
these discrepancies when the LLL is taken into account. 

Following Reikuba (2001), we have assumed a distance 
modulus of 27.8 and we have corrected for extinction the ob- 
served photometric bands with a mean E(B - V) = 0.1 us- 
ing A v - 4.40, A v = 3.1, A K = 0.38, according to th e val- 
ues quoted in the ADPS project (Moro & Mun arifcOOoT) 11 . We 
have grouped the GCs according to their L c £' s luminosities: 
L dusi L LLL < 5 ^ L cIu S/l lll E [5i 10]i L dus fjiXL 6 [io,30], and 

L dus/ L LLL > 3 We have use( j the L LLL yame at j Q yJ . ^ 

corresponds to M™" = -8.5 mag. The resulting color-color 
diag ram i s show n in Fig. |6] The results of SSP models from 
Worthed ill 994 have also been plotted for comparison. Note 
that a rigorous comparison should include the confidence inter- 
vals for the results of SSP models; however, this would require 
knowledge of the correlation coefficients between the different 
luminosities, which, unfortunately, have onl y been evaluated 
for the case of young stellar populations (see ICervinoll2003dt 
ICervifio et al.l2002bl) . 



11 Available at http : //ulisse . pd. astro . it/Astro/ADPS/ 
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The figure shows several interesting results. First of all, the 
cluster sample tends to have V-K colors redder than the pre- 
dictions of SSP models, as is expected from a luminosity lim- 
ited survey. Second, there is a statistical trend of the position in 
the plane with the Lk value: (i) The first bin of clusters, the one 
with L°£ S /L^ LL < 5, contains five of the 23 clusters plotted, and 
one of them is the cluster that disagree the most with the SSP 
results, (ii) Of the twelve clusters with Lf' s /L^ LL e [5, 10], 
five disagree with the results of synthesis models, (iii) There 
are four clusters with L C ^ US /L^ L e [10,30], and only one of 
them falls far from the results of SSP models, (iv) The two 
clusters with U^' s IL L ^ L > 30 are reproduced reasonably well 
by SSP models. In summary, clusters with lower Lk tend to de- 
viate more from SSP models than luminous clusters, as can be 
expected when sampling effects are present. 

5. Caveats about self-consistency and limitations 

In the present work we have established a luminosity limit, 
the LLL, for the application of synthesis models to observed 
clusters. The definition we chose for the LLL might raise ob- 
jections, in at least two different ways. First, the LLL is im- 
plicitly defined by the isochrones and the atmosphere libraries 
assumed in the code; a natural doubt is therefore whether as- 
suming different input libraries in the code would change our 
results. Second, there are evolutionary phases that are char- 
acterized by extremely high luminosities and short lifetimes. 
According to our definition, it is these phases that determine 
the value of the LLL, at the ages they occur; nevertheless, they 
have a very small probability of being observed in real clusters. 
Therefore we face the apparent oddity of setting a limit much 
more restrictive than necessary. 

This section is devoted to clarify these points and to con- 
vince the reader that our definition of the LLL is operationally 
good and, to a satisfactory extent, physically based. We will 
also discuss its residual degree of arbitrariness, and the limi- 
tations it implies. Finally, we will also explain why, by virtue 
of its limited reach, the exact definition of the LLL is not very 
crucial in practice. 

5.1. Dependence on input isochrones and model 
atmospheres 

A possible objection to our definition of the LLL is that 
it depends on the input libraries (isochrones and atmo- 
sphere models) assumed by the synthesis code, suggesting 
that it is not physically grounded. For example, the LLL can 
be lowered if some stages of stellar evolution are not in- 
cluded in the computations of the synthesis code, as is of- 
ten the case for the Asymptotic Giant Branch (AGB) phase, 
which is ne glected, e.g., in starbu r st-orie nted codes such as 
the code bv ICervifio & Mas-Hessd d 19941) and starburst99 12 
jLeitherer et alJll999l) . Obviously, the results of evolutionary 
synthesis models become intrinsically incomplete at the ages 
where the evolutionary stages neglected are relevant; in the case 
of the AGB phase, ages older than 50 Myr. Therefore, under 

12 Available at http : //www . stsci . edu/science/starburst99/ 
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these circumstances the concept of LLL cannot be applied, not 
because it is ill-defined, but rather because the code should not 
be used in the first place. 

Furthermore, the agreement shown in Fig. @] between 
the trend of our predictions for M' mn and the Atio% values 
by various authors, which have been obtained with different 
isochrones and atmosphere libraries, suggests that the LLL 
value does not substantially depend on the input libraries as- 
sumed - provided they include the same evolutionary stages. 
This conclusion should not surprise too much, since the results 
of different theoretical models in the fields of stellar evolution 
and atmosphere modeling are substantially converging. It must 
be noted, however, that it should not be of great concern even 
if this were not the case, since at the LLL level one must only 
worry about self-consistency, i.e. that the LLL is taken from the 
same input libraries that enter the synthesis code; the judgment 
on the input data quality has already been made with the choice 
of the synthesis code. Stated otherwise, the choice of a synthe- 
sis model is logically prior to the computation of the LLL, and 
it already implies trusting the input data the model is based on. 

5.2. Dependence on short, very luminous evolutionary 
phases 

Some readers may have noted that our definition does not 
involve stellar lifetimes, but just their luminosities. Such a 
definition might seem too restrictive in those evolutionary 
phases characterized by very luminous but short-lived ob- 
jects: Supernovae (SNe), luminosity peaks of Luminous Blue 
Variable stars, and recurrent He-shell flashes in TP-AGB stars 
are cases in point. Due to their short duration, the probability of 
observing such phases in real clusters is vanishingly small, and 
they do not usually contribute to the observed luminosities. Yet, 
from a theoretical standpoint these phases provide the maxi- 
mum luminosity: hence they determine the LLL value at the 
corresponding ages. Therefore, it may seem that our definition 
is at risk of giving too much weight to luminous, short-lived 
evolutionary phases; if this were the case, the definition would 
be too restrictive, and it would mistakenly ban the application 
of synthesis models to cases where they would otherwise be 
applicable. Although correct in principle, this conclusion does 
not hold in practice, as we will make clear in the following. 

First, very short evolutionary phases are not usually in- 
cluded in isochrones, due to the same reason for which they 
should not determine the LLL: they are not expected to be 
observed. Of course, this argument would no longer be valid 
should these rapid phases be included in isochrones: but users 
of synthesis models should always take the elementary caution 
of documenting themselves about the properties, assumptions, 
and limitations of the models they use. 

Second, it is not easy finding cutoff criteria alternative to 
the LLL and equally easy to apply, to discriminate systems with 
poor statistics from systems with good or unknown statistics. 
An alternative definition could be based, for example, on the 
comparison of the cluster luminosity with the sum of the lu- 
minosities of all the individual stars included in the isochrone 
(rather than the most luminous one alone, as in our defini- 



tion). Unfortunately, such criterion depends on the mass res- 
olution - the number of points - of the isochrone, which in turn 
depends on the individual, idiosyncratic interpolation philoso- 
phy of each code. Therefore, such definition would yield pro- 
foundly different results for different codes, even if the same 
set of isochrones and atmosphere libraries were adopted. Still 
another definition could rely on the energy emitted by a star 
during a particular phase, rather than its luminosity, in order to 
avoid the problem that arises with very luminous but short stel- 
lar phases. However, it can be seen that in this case the limit 
depends on how evolutionary phases are defined. For example, 
two possible alternatives would be either (i) "standard" evolu- 
tionary phases (e.g., Wolf-Rayet, RSG, AGB stars, and so on), 
or (ii) the interval around each point on the isochrone. In ei- 
ther case, the numerical result would depend on the adopted 
definition of evolutionary phase, hence on an arbitrary choice. 

As a final remark, it should be emphasized that neither the 
LLL nor Af™ define the statistics and the confidence levels 
of synthesis models, which are defined, for example, by the N 
formalism. In this sense, the LLL is just a two-state switch that 
indicates whether in each given case a statistical approach is 
mandatory; if it is not, it is still possible that sampling effects 
do play a role. In epidemiological terms, the LLL criterion is 
extremely specific, but not very sensitive: it does not give false 
positives ("good" clusters with rich statistics misunderstood for 
"bad" clusters with poor statistics), but it might give many false 
negatives ("bad" clusters with poor statistics not identified). As 
such, it gives a much more basic information than N, hence it 
is not too worrisome whether its definition has a certain degree 
of arbitrariness. Ideally, the predictions performed by any syn- 
thesis code should be accompanied by a statistical analysis of 
the system under study, bypassing the need to resort to coarse 
indicators like the LLL. 

6. Summary and conclusions 

In this paper, we have drawn the attention to a basic condition 
that must be fulfilled in order to interpret observed stellar clus- 
ters by means of synthesis models: 

The total luminosity of the cluster modeled must be larger 
than the individual contribution of any of the stars included in 
the model. 

In particular, the luminosity of the most luminous star in- 
cluded in the model defines a Lowest Luminosity Limit below 
which a cluster suffers from severe sampling effects, and it can- 
not be modeled by an analytical synthesis code. We have em- 
phasized that, although our method does not recognize all the 
clusters affected by sampling effects, by virtue of its simplic- 
ity it can be used to separate out the most critical cases from 
large surveys without the need to perform a sophisticated and 
time-consuming statistical analysis. Transient events must be 
considered in the LLL only if they are included in the results of 
SSP computations; if they are not included, the definition of the 
LLL is fully meaningful, under the assumption that such events 
are not present in the observations. Although this assumption 
might seem, at first glance, a very restrictive one, it does not 



10 



M. Cervino and V. Luridiana: The Lowest Luminosity Limit 



in fact limit the application of synthesis models, since the SSP 
results themselves are also only valid under the same assump- 
tion. 

We have obtained the Lowest Luminosity Limit for a wide 
range of ages and metallicities. In turn, the Lowest Luminosity 
Limit has been used to obtain the associated mass limit A\ mn , 
which is the cluster mass below which the colors may be 
severely biased with respect to the results of synthesis models. 
The limit depends on the age and the metallicity of the cluster 
and is more stringent in the near infrared bands; in particular, 
it has a value between 10 4 and 10 5 M G for K. We have also 
shown that the luminosities in different bands of clusters with 
mass ten times larger than M mm have relative dispersions of 
10% at least: that is, yVf""" can be used to estimate a lower limit 
to the expected dispersion of luminosities 

As observational applications, we have shown that in sur- 
veys that reach luminosities near the Lowest Luminosity Limit 
the color distributions will be skewed toward the luminosity 
with the lowest N ssp value, and that the skewness may be used 
to constrain the distribution of initial cluster masses M, pro- 
vided the other relevant cluster parameters are known. We have 
also analyzed a sample of Globular Clusters in NGC 5128, 
showing that sampling effects are probably playing a role, and 
that they can explain, at least partially, the bias of the observed 
colors with respect to the predictions of synthesis models. 

Finally, we have discussed the assumptions underlying our 
definition of the LLL, its virtues, and its drawbacks. In particu- 
lar, on the positive side we have stressed the extreme simplicity 
of the LLL as a criterion to recognize in an economical way the 
occurrence of poor statistics in observed clusters; on the nega- 
tive side, we have highlighted the limited capability of the LLL 
criterion to recognize all the cases of poor statistics. 

The preceding results have been obtained from very ba- 
sic concepts implicit in the modeling performed by synthesis 
codes. In particular, the definition of the LLL is very simple, 
and its computation can be as elementary as finding a maxi- 
mum in a table, at least in the case of SSP models. But this 
simplicity may be deceiving, for simple as the concept may be, 
it has far-reaching consequences, which have been too often 
overlooked in the analysis of stellar clusters. Furthermore, al- 
though the simplicity of the concept presented here might dis- 
appoint those who are fond of spectacular or complicated re- 
sults, we consider it as an asset rather than a drawback, since 
the simplicity of our method may serve as an incentive to those 
who realize the importance of sampling effects but do not know 
how to deal with them. Therefore, our aim in the present work 
is twofold and goes beyond the particular results we have de- 
scribed in this paper: in the short term, we want to provide a 
very simple sieve to detect rapidly clusters with poor statis- 
tics; in the long term, we hope to motivate other researchers, 
both theoreticians and observers, to move in this direction and 
explore the implications of sampling effects in the analysis of 
the integrated light from clusters. In particular, we urge all the 
model makers to include a proper statistical treatment in their 
code, and all the model users to take into account the effects 
of sampling in the interpretation of the data. Current synthesis 
models are an optimum tool for the interpretation of the aver- 



age properties of stellar systems, but they need to be updated to 
keep pace with the new observational data. 
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Fig. 1. LLL in magnitudes for the bolometric luminosity and different bands, as a function of the age for three different metal- 
licities: top panels Z=0.019 (solar); middle panels Z=0.008; and bottom panels Z=0.004. Left panels correspond to isochrones 
with a simplified TP-AGB evolution, and right panels to isochrones with a detailed TP-AGB evolution. The different lines corre- 
spond to the following: bold solid line: bolometric magnitude; solid line with triangles: U; dashed line: B; dot-dashed line with 
triangles: V; dotted line: R; dash-dot-dot-dotted line: /; solid line: J; dashed line with triangles: H; and dot-dashed line: K. 
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Fig. 2. LLL in magnitudes for the bolometric luminosity and different bands, as a function of the age for Z=0.001 (top left), 
Z=0.0004 (top right), Z=0.0001 (middle left) and Z=0 (middle right). Symbols like in Fig.^ The bottom panel compares the 
LLL for V and K for Z=0.019 and Z=0.0001 . 
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Fig. 3. M mm for the bolometric luminosity and different bands, as a function of the age and the metallicity. Symbols as in Fig.^ 



